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ABSTRACT 


The  least  constant  density  of  a  spherical  swarm  of  particles  is  first  reviewed 
for  the  case  of  circular  motion  about  a  fixed  attracting  mass.  The  density  ratio 
is  defined  as  the  ratio  of  the  density  of  the  swarm  to  the  attracting  mass  divided 
by  the  volume  of  a  sphere  whose  radius  is  the  radius  of  the  circular  orbit.  The 
classical  linearized  analysis  shows  that  the  density  ratio  must  be  at  least 
three.  This  work  is  extended  to  take  into  account  the  finite  size  and  the  variable 
density  of  the  swarm;  the  least  density  ratio  is  increased  in  each  case.  The 
problem  is  then  generalized  to  the  case  where  the  swarm  is  located  at  either 
of  the  stable  Lagranglan  equilateral  triangle  points  (L^,Lg)  in  the  restricted 
problem  of  three  bodies.  The  least  density  ratio  is  now  3  -  (9/4)|ji  to  first  order 
terms  in  (jl,  and  the  mass  is  that  of  the  combined  earth-moon  system;  since  p 
is  approximately  1/82  in  the  earth-moon  system,  the  change  in  density  from  the 
classical  value  is  slight. 

The  recent  (1961)  astronomical  observations  of  K.  Kordylewski  of  Krakow 
Observatory  (Poland)  with  respect  to  the  existence  of  a  pair  of  cloud-like 
satellites  near  each  stable  L-point  are  discussed,  and  a  dynamical  model  is 
analyzed  which  yields  results  that  are  in  substantial  agreement  with  the 
observations.  The  model  assumes  the  existence  of  a  small  mass  at  either  stable 
L-point  and  dust  particles  moving  in  the  surrounding  region.  The  main  tool  in 
the  analysis  is  the  Jacobi  (energy)  integral  which  determines  the  surfaces  of 
zero  relative-velocity;  the  model  indicates  that  for  an  appropriate  and  plausible 
range  of  parameters  the  surfaces  of  zero  relative-velocity  can  split  into  shapes 
that  are  similar  to  those  noted  by  Kordylewski. 
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INTRODUCTION 

Studies  in  the  nebular  hypothesis  and  in  the  theory  of  comets  have  led  to  the 
consideration  of  the  stability  of  a  spherical  swarm  of  particles  in  circular 
motion  about  a  massive  center  of  attraction  (sun).  Any  particle  of  the  swarm 
is  attracted  toward  the  sun,  and  toward  the  center  of  the  swarm  by  particles 
interior  to  it.  If  it  is  assumed  that  the  particles  are  at  least  100  microns  in 
diameter  [1],  we  can  restrict  ourselves  to  gravitational  forces  alone  and  ignore 
radiation  forces  as  well  as  electrical  ones. 

To  illustrate  the  dynamical  situation  in  the  simplest  way  consider  a  pair 
of  equal  small  particles  (m^sm^^m)  in  circular  motion  about  a  fixed  center  of 
attraction  M,  (m/M  <<  1),  as  shown  in  Figure  (1).  For  simplicity  we  assume 
that  the  orbits  of  m^  and  m^  are  in  the  same  plane.  The  particles  m^  and  m^ 

can  move  in  circular  orbits  of  radii  (b-a) 

and  (b+a),  if  the  sole  fcxce  on  either 

is  due  to  M  .  The  respective  angular 

velocities  of  the  radii  vectors  are  then 
2  3  2 

given  by  =  GM/ (b-a)  and  =  GM/(b+a) 

Figure  1 

Sponsored  by  the  Mathematics  Research  Center,  United  States  Army,  Madison, 
Wisconsin  under  Contract  No.:  DA-11-022-ORD-2059. 
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where  G  is  the  gravitational  constant.  Thus  ^  and  in  fact  • 

If  the  particles  are  in  the  same  radius  at  time  t  =  0,  the  inner  particle  will 
move  more  rapidly  than  the  outer  one  and  the  particles  will  separate;  the 
distance  between  them  will  be  >  2a  and  Increasing  until  they  are  180°  apart. 

However  if  the  forces  due  to  mutual  attraction  between  m^  and  m^  are 
included  in  the  analysis,  the  force  on  m^  will  be  diminished  and  that  on  m^ 
will  be  increased  at  t  =  0  .  This  action  will  decrease  and  increase 
and  it  is  possible  to  make  proper  choice  of  the  parameters.  It  is 

only  necessary  to  require  that 

GMm/(b-a)^  -  Gm^/(4a^)  =  mw^(b-a), 

GMm/(b+a)^  +  Gm^/  {4a^)  s  mw^{b+a), 

if  the  particles  are  to  remain  on  the  same  radius  to  M  for  all  time.  If  we 
introduce  a  dimensionless  ratio  of  frequent  use,  to  be  called  the  density 
ratio  X  ,  as 

.  .  f  (Zmlirb^  1 .  r__2m__  IfUlZabil 

■  lua)’Jl“J'L(4./JXZa^  Jl  “  J 


then  elimination  of  w  between  the  two  equations  yields 


_  bisbii^ 

"(b^-aV 


(1) 


The  quantity  K  is  a  ratio  of  densities.  The  first  term  is  the  total  particle  mass 
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(2m)  divided  by  the  volume  of  a  sphere  of  radius  (2a);  the  second  term  can 
be  interpreted  as  the  mass  of  the  primary  (M)  divided  by  the  volume  of  the 
orbit  (radius  b).  If  b»  a,  which  means  that  the  distance  between  the 
particles  is  small  compared  with  b,  then  X.  =  3  very  nearly.  The  density 
ratio  must  therefore  be  3  for  the  assumed  motion  to  persist.  If  the  value  of 
X  is  <  3  then  the  self  attractive  forces  are  insufficient  to  keep  the  particles 
together  and  the  simple  model  shows  that  the  density  ratio  of  a  swarm  of 
particles  must  exceed  a  minimum  value  for  the  swarm  to  remain  together. 

The  simple  two-particle  model  is  not  realistic  enough  for  the  discussion 
of  a  swarm  of  particles  more  nearly  spherical  in  distribution.  In  fact  the 
assumed  motion  is  unstable;  the  application  of  this  result  is  of  some  interest 
in  other  problems  (orbital  rendezvous  and  docking  dynamics  of  space  vehicles) 
and  is  considered  further  in  the  Appendix. 

The  self-attractive  forces  in  a  spherical  swarm  are  of  a  more  complex 
nature  than  in  the  two-particle  model;  the  force  exerted  on  any  particle  of  a 
spherical  swarm  with  radially  symmetric  density,  by  the  remainder  of  the 
swarm,  will  tend  to  zero  as  the  particle  approaches  the  center  of  the  swarm. 

The  self  attractive  force  can  become  arbitrarily  large  in  the  two-particle  model 
and  this  singularity  makes  the  model  unsuitable  for  further  study  of  the 
stability  of  extended  bodies  composed  of  small  particles. 

The  classical  analysis  of  a  spherically  symmetric  swarm  of  particles  assumes 
that  the  radius  of  the  swarm  is  infinitesimally  small  compared  with  the  radius 
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of  the  orbit  and  that  the  density  is  constant.  This  work  is  briefly  reviewed  in 
I  and  extended  to  the  case  of  a  finite  swarm  radius  by  means  of  an  energy 
integral.  The  influence  of  variable  density  can  also  be  taken  into  account 
in  an  approximate  way  under  reasonable  simplifications.  Both  effects  Increase 
the  least  density  ratio. 

It  is  of  Interest  to  extend  the  basic  result  of  the  linearized  analysis  to  the 
case  of  motion  with  respect  to  two  finite  bodies  in  the  sense  of  the  restricted 
problem  of  three  bodies.  In  this  formulation  the  two  finite  bodies  rotate 
unifcarmly  in  circles  about  their  common  mass-center  and  the  third  body  is 
infinitesimal  in  mass.  That  is,  the  third  body  moves  under  the  attractive 
forces  of  the  two  finite  bodies,  but  does  not  influence  their  motion.  The 
spherical  swarm  of  particles  is  the  third  (infinitesimal)  body. 

It  is  well  known,  in  the  theory  of  the  restricted  problem  of  three  bodies, 
that  five  equilibrium  positions  are  possible,  [2,  chap.  8]  .  The  three  straight- 
line  configurations  (denoted  by  Lj,  L^,  L^)  are  unstable  for  any  mass-ratio 
of  the  finite  bodies.  The  two  equilateral  triangle  solutions  (denoted  by  L 

4^ 

Lg)  are  unstable  if  ^  >  0, 0385  where  n  is  the  raUo  of  the  smaller  finite 
mass  to  the  total  mass.  The  L^,  Lg  points  are  the  vertices  of  equilateral 
triangles  with  the  two  finite  masses  at  the  remaining  vertices  and  in  the  plane 
of  the  orbit.  Since  a  small  parUcle  placed  at  L^,  Lg  in  the  earth-moon  system 
(jx  <  0. 0385}  is  at  a  posiUon  of  stable  equilibrium,  at  least  according  to  the 
small  displacement  criterion  of  a  linear  analysis,  it  will  tend  to  remain  in  this 
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vicinity  for  a  long  time;  ultimately  the  particle  will  leave  the  L^,  region, 
to  be  replaced  by  other  particles  that  are  moving  in  the  surrounding  space.  It 
is  therefore  possible  to  have  an  accumulation  of  matter  at  the  L^,  Lg  points. 

However,  the  condensation  of  such  matter  into  a  stable  swarm  would 
again  require  considerations  similar  to  that  in  I  for  a  swarm  moving  in  a 
circular  orbit  about  a  fixed  mass.  In  the  present  case  the  motion  of  a  swarm 
centered  at  or  Lg  would  again  be  circular  about  the  mass -center  of  the  two 
finite  bodies.  The  precise  considerations  for  this  case  are  given  in  II  for 
the  case  of  a  spherical  swarm  of  constant  density. 

In  the  analyses  outlined  so  far,  the  typical  particle  of  the  swarm  was 
assumed  to  be  internal  to  it.  This  is  necessary  for  the  discussion  of  the 
stability  of  the  swarm  and  the  determination  of  the  least  value  of  the  density 
ratio  \  .  However,  it  is  also  necessary  to  consider  a  particle  external  to 
the  swarm.  This  is  the  case  in  problems  which  involve  the  tall  of  a  comet. 

Here  we  are  concerned  with  particles  external  to  the  nucleus  of  a  comet;  the 
tail  is  assumed  to  be  so  tenuous  that  the  particles  in  it  do  not  attract  each 
other  but  are  influenced  solely  by  the  nucleus  and  the  masses  which  can 
influence  the  motion  of  the  nucleus. 

These  considerations  bring  us  to  the  problem  treated  in  III,  in  relation 

to  the  remarkable  observations  made  at  the  Krakow  Observatory  by  the  Polish 

astronomer  K.  Kordylewski.  In  the  search  for  additional  natural  satellites  of 

the  earth,  Kordylewski  looked  for  such  objects  near  L.  and  L-  .  His 

4  5 
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observations  indicate  that  two  thin  cloud^like  satellites  exist  at  each  point. 

The  problem  treated  in  III  is  this:  is  it  possible  to  explain  these  observations 
in  terms  of  a  simple  dynamical  model  ? 

The  model  assumes  the  existence  of  a  small  nucleus  at  either  L  point 
(L^  or  Lg)  and  the  existence  of  tenuous  matter  near  it.  As  a  simplifying 
assumption  the  action  of  the  moon  is  neglected  and  the  nucleus  is  assumed  to 
move  in  a  circular  orbit  about  the  earth  at  a  distance  equal  to  that  of  the  moon 
from  the  earth.  It  is  a  simple  matter  to  derive  a  Jacobi  integral  for  the 
particles  external  to  the  nucleus.  If  the  velocity  of  the  particles  is  very  close 
to  the  orbital  speed  at  that  point  (L),  it  can  be  shown  that  the  surfaces  of 
zero  relative  velocity  become  two  small  elongated  tube-like  surfaces  near  each 
L  point.  If  the  relative  velocity  increases,  these  surfaces  are  much  enlarged 
physically  or  cease  altogether,  and  the  particles  cannot  remain  near  the 
nucleus  for  any  length  of  time.  Under  appropriate  conditions  not  only  can  the 
particles  remain  near  the  zero  relative  velocity  surfaces,  but  if  they  do,  they 
must  do  so  near  two  small  regions.  Kcx'dylewski's  observations  at  each  L 
point  have  Indicated  just  such  cloud-like  formations;  the  observational 
problems  are  serious  ones  but  observations  [14,  15,  16,  17]  have  been  announced 
fcr  each  L-point.  Considerable  work  has  also  been  done  by  others  in  the 
search  for  natural  satellites  which  are  much  nearer  to  the  earth.  Here  again 
the  observaUonal  difficulties  are  great  but  of  a  different  type;  no  announcements 
have  been  made  as  to  the  discovery  of  such  near-earth  natural  satellites. 
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LIST  OF  MAIN  SYMBOLS 

distance:  orbital  radius;  semi-major  axes 
eccentricity  of  orbit  {positive) 

3  -1  -2 

gravitational  constant  (units:  L  M  T 
constant  (twice  rate  of  change  of  area  swept  out) 
constant  In  Jacobi  Integral 
(units  of  k:  T 

Lagrangian  equilibrium  point,  in  particular  either  L  or  L_ 

4  5 

Masses 

distances 

coordinates 

coordinates 

relative  velocity  in  Jacobi  integral 
true  anomaly 

3 

constant  (mass/radius  of  swarm) 

(3^1/4)  0-2^1) 

ratio  of  smaller  finite  body  to  the  sum  of  both  finite  bodies  in 
the  restricted  problem  of  three  bodies,  (0<  |x  ^  1/2) 
density  ratio,  or  the  ratio  of  the  density  of  the  swarm  to  the 
mass  of  the  primary  divided  by  the  volume  of  a  sphere  whose 
radius  is  the  radius  of  the  orbit, 
coordinates 

kt,  dimensionless  time 
angle 

polar  coordinates 

angular  speed 

d/dt  (time  differentiation) 

d/dr  (dimensionless  time  differentiation) 


The  subject  of  interplanataxy  dust  has  been  considered  in  several  recent 
(and  many  past)  studies,  [1,  3,  4]  .  The  present  requirements  are  reasonable 
ones  since  such  small  dust-like  particles  are  continuously  moving  in  the 
spaces  of  the  solar  system.  In  fact  the  discussion  of  the  phenomenon  of  the 
Gegenschein  [2,  p.  305]  and  [9]  has  some  points  of  similarity  to  the  present 
case  even  though  the  L^,  points  are  stable  positions  of  equilibrium;  in  the 
Gegenschein  problem  the  L  points  (L^,  L^)  are  unstable  positions  of 

equilibrium.  However,  Kordylewski' s  observations  indicate  that  a  pair  of 
cloud-like  satellites  do  exist  near  each  stable  L-point.  If  it  is  assumed  that 
some  matter  can  accumulate  at  L.  or  L.  in  the  form  of  a  nucleus,  the 
dynamical  model  can  explain  the  form  of  the  clouds  observed  by  Koardylewski. 
Finally,  the  least  density  of  the  nucleus,  considered  as  a  swarm,  is  of 
possible  Interest  in  connection  with  space  probes  near  the  L.,  L.  points;  the 
actual  density  would  seem  to  be  comparable  with  the  density  of  the  atmosphere 
of  the  earth  well  inside  the  re-entry  altitudes  (less  than  50  miles). 
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I.  CIRCULAR  MOTION  ABOUT  A  FIXED  MASS-POINT. 


(a)  Basic  Equations;  Linearized  Analysis. 

In  Figure  (2)  consider  the  mass  M  to  be  fixed  at  the  origin  0  of  the 

(x,  y,  z)  coordinate  system.  The  center  of  the  spherical  swarm  is  at 

C(|,  T^)  in  the  z  =  0  plane  and  point  C  moves  uniformly  in  a  circle  of  radius 

2  2  2 

a ,  where  a  =  |  +  t)  ,  and  r  is  the 

distance  of  any  point  P(x,  y,  z)  of  the 

swarm  to  the  origin  0,  where 
2  2  2  2 

r  =  X  +  y  +  z  .  The  swarm  is  assumed 
to  be  spherical  in  shape  and  constant  in 
density.  If  or  is  the  constant  density 
of  the  swarm  then  the  mass  contained  in 

3 

a  sphere  of  radius  PC  is  m  =  (4ir/3XPC)  o- 

3 

or  m/(PC)  =(4/3)Tro-  =  p  =  constant. 

Under  these  assumptions  it  is  easy  to  write  the  differential  equations  of  motion 
of  P  with  respect  to  the  inertial  system  (x,  y,  z)  as 


X  =  -GMx/r^  -  GP(x  -  |)  , 

y  =  -GMy/r^- Gp{y  -  Ti)  ,  (2) 

z  =  -GMz/r^  -  Gpz  , 


where  ,  =  d/dt  and  G  is  the  gravitational  constant. 
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The  general  equations  (2)  hold  at  P  =  C,  that  is,  at  x  =  |,  y  =  ti, 
z  s  0  and  r  =  a  =  constant: 

i  =  -GM|/a^,  1^  =  -GMri/a^  .  (3) 

2  3  "2 

If  k  =  GM/a  the  equations  (3)  both  assume  the  simple  form  |  +  k  §  =  0 

and  we  choose 


I  =  acos(kt)  , 

T|  =  asin(kt)  , 


(4) 


as  the  coordinates  of  C  as  a  function  of  time  t  .  The  equations  (4)  show 
that  it  is  convenient  to  introduce  kt  s  t  as  the  dimensionless  time  and  a  as 
the  new  unit  of  length;  the  dot  ( ' )  above  a  quantity  denotes  differentiation  with 
respect  to  the  time  t  and  a  prime  ('}  above  a  quantity  denotes  differentiation 
with  respect  to  the  dimensionless  time  kt  =  T:k(*}=(  }'  .  We  also 
introduce  new  variables  for  (x,  y,  z}  in  a  rotating  reference  system  in  which 
C  is  the  fixed  point  at  the  x-axis  at  unit  distance  hrom  0  : 

X  =  a[(u-t'l}cos  T  -  vsinr] 

y  =  a[  (u4-l)sin  x  +  vcosr]  (5) 

z  =  aw  . 

Thus  u,  V,  and  w  are  all  small  quantities  compared  with  unity.  If  the 
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substitutions  (4)  and  (5)  are  made  in  Eqs.  (2)  and  the  trigonometric  terms 
are  eliminated  the  results  can  be  written  as 

u"  -  2v'  -  u  =  1  -  (u+l)/R^  -  Xu  , 

v"  +  2u'  -  V  =  -v/R^  -  Xv  ,  (6) 

3 

w"  =  -w/r  -  Xw  , 

2  3  3 

if  use  is  made  of  k  =  GM/a  ,  and  X  =  [(4u'/3)a  /M]<r  is  the  appropriate 

2  2  2  2 

density  ratio,  and  R  =  1  +  2u  +  (u  +v  +w  )  .  In  the  linearized  analysis, 

such  as  that  outlined  in  Routh  [6,  p.  264-267]  it  is  assumed  that  R*"^  s  l-2u 

>3 

and  R  =  l-3u  and  the  linearized  equations  are  then  obtained  in  the  form 
from  Eqs.  (6)  as 

u"  -  2v'  +  u(X  -3)  =  0  , 


v"  +  2u'  +  v(X)  =  0  ,  (7) 

w"  +  w(l  +  X )  =  0  , 


if  squares  of  u,  v  or  w  are  neglected  in  the  expansions.  The  third  of 
equations  (7)  has  only  oscillatory  solutions  since  the  density  ratio  is 
positive.  If,  in  the  first  two  equations  we  set  u  =  Ue’'’^  and  v  *  the 

condition  on  v  is 

|[v^  +  (X-3)]  (-2v)  I 


(2v) 


[v^+X] 


=  0  , 


(8) 
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or 

+  (2X.  +  l)v^  +  (X^  -  3K)  =  0  . 

2 

For  stability  v  must  be  negative;  since 

2v^  =  -{2X  +  1)  A  n/i  +  16X  , 

2 

the  requirement  is  that  (1  4- 16X)  be  less  than  (1  +  2X}  since  X  is  positive. 
This  leads  to  the  known  condition 

X  >  3  (9) 

which  is  the  classical  one  discussed  in  Routh  [6,  p.  266]. 

According  to  Tisserand  [  5,  p.  258£f],  Schiaparelli  obtained  the  criterion 
on  statical  grounds  that  if  X  <  2  the  swarm  would  dissolve  due  to  the  action 
of  the  attracting  mass.  Routh's  statement  [  6,  p.  264]  that  Schiaparelli 
required  X  >2  for  stability  is  not  quite  accurate  since  only  the  weaker,  but 
correct,  assertion  was  made  that  the  swarm  would  dissolve  for  X  <  2  .  The 
dynamical  argument  shows  that  the  swarm  would  also  dissolve  if  X  <  3  . 

Tisserand  [  5,  p.  269ff]  also  discussed  the  dynamical  arguments,  which 
are  due  to  Charller  and  Picart;  these  are  presented  in  simplified  form  in 
Routh  [6,  p.  264f].  The  case  where  the  orbit  is  non-circular  was  also  noted 
briefly  by  Tisserand  [5,  p.  275]  in  connection  with  Floquet's  theory  of 
ordinary  differential  equations  with  periodic  coefficients.  Tisserand  noted  that 
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if  the  eccentricity  e  of  the  orbit  is  small  compared  with  unity,  then  \  >  3  +  He 
but  he  did  not  evaluate  the  constant  H  .  Routh  [6,p.  266]  has  given  an 
approximate  treatment  to  show  that  H  =  5  ,  The  non-circularity  of  the  orbit 
therefore  raises  the  least  value  of  the  density  ratio.  Routh,  [6,p.  406]  has 
also  discussed  the  equally  Important  case  where  the  swarm  is  non-spherical 
in  shape  but  the  orbit  is  circular.  Picart  [  7]  emphasized  the  importance  of 
using  the  true  anomaly  as  the  new  independent  variable  in  those  cases  where 
the  orbit  was  non-circular;  see  also  [10].  Callandreau  [8,  11,  12]  also 
contributed  to  the  study  of  this  question.  Some  Indication  of  the  value  of  the 
transformations  used  by  Picart  will  be  noted  briefly  in  lie  . 

It  is  also  possible  to  approach  the  problem  of  non-circular  motion  in  a 
different  manner  although  we  shall  not  go  into  details.  For  example,  we  can 
assume  that  a  spherical  swarm  is  in  aii>itrary  elliptical  (parabolic,  hyperbolic) 
motion  about  a  fixed  M  ;  consider  a  particle  on  the  surface  of  the  swarm  and 
acted  upon  by  the  particles  internal  to  it,  by  M,  and  by  assumed  surface 
forces  (tangential  and  normal)  which  constrain  the  particle  to  remain  on  the 
surface.  Thus  the  motion  of  the  particle  is  given  and  the  equations  of  motion 
determine  the  value  of  the  surface  forces.  The  normal  surface  force  can  now 
be  set  -equal  to  zero.  In  the  case  of  a  circular  orbit  this  condition  again  leads 
to  the  requirement  \  =  3  and  this  value  increases  with  e  . 
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(b)  Tacobl  (Eneravl  Intecrral;  Arbltraarv  Swarm  Radius. 

The  governing  non-linear  differential  equations  [6]  can  be  integrated 
once  if  the  equations  are  multiplied  respectively  by  u* ,  v' ,  w'  to  get 

(u')^  +  (v')^  +  (u^+v^)(X  -  1)  -  2u  -  "■  ■"  ^  ■■  —  +  2  =  const  (10) 

'vtl+u)^+  v^ 

where  the  constants  have  been  appropriately  adjusted  and  the  dependence  upon 
w  has  been  assumed  to  be  negligible  for  simplicity.  The  constant  is  zero 
if  u  =  V  =  0  . 

If  u  =  rcos(0),  V  =  rsin{e)  and  if  we  write 

(u*)^  +  (v')^  ®(r,  0;  X)  s  const.  (11) 

then  «  (r,  0;  X)  =  (X  -  l)r^  -  2rcos(0)  -  .=.-=:rL  ;  +  2  ,  (12) 

VI  +  2rcos0  +  r 

The  least  value  of  X  is  required  such  that  #  >  0  for  given  (u,  v)  or  (r,  0)  , 

where  (0  <  r  <  1)  .  For  example  if  u,  v  are  small  and  only  quadratic  terms 

2  2 

are  retained,  «  =  (X  -  3)u  +  v  and  «  >  0  if  and  only  if  X  >  3  for 

aibitrary  u,  v;  this  is  the  known  linear  criterion. 

In  the  non-linear  case  we  can  first  assume  that  r  is  a  fixed  quantity  and 

seek  to  extremize  4;  this  fixes  0  as  a  function  of  r  and  three  cases  can 

arise;  cos  0  =  -r/2  or  sin0  =  0  (cos  0  =  +1  or  -1)  .  The  first  case  leads  to 
2 

«  =  X  r  which  does  not  restrict  X  .  In  the  two  following  cases  the  condition 


cos  6  =  -1  leads  to  the  larger  value  of  \  in  the  form 


X,  ,3tit3z2rl 

tl-i) 

and  this  condition  dominates  the  problem  and  *  attains  its  minimum  value  under 
these  conditions.  If  r  is  assumed  to  be  negligibly  small  with  respect  to 
unity  then  the  linear  result  X.  =  3  is  again  obtained.  Otherwise  X  =  3+3r  to 
lowest  order  terms  in  r;  in  general  the  density  ratio  must  increase  if  the 
swarm  increases  in  size.  This  is  a  result  to  be  expected  on  physical  grounds 
and  indeed  we  would  expect  the  result  to  become  unbounded  as  r  =  1  . 

(c)  Swarm  of  Variable  Density. 

Since  the  radius  of  the  swarm  is  almost  always  small  compared  with  the 
radius  of  the  orbit,  the  effect  of  the  finite  size  on  X  is  rather  small.  However, 
the  effect  of  variable  density  is  likely  to  be  more  marked.  To  take  approximate 
account  of  this  effect  we  return  to  the  first  two  non-linear  equations  (6)  and 
seek  to  modify  the  equations  by  a  reasonable  assumption  about  the  variation 
of  density  with  the  radius  of  the  swarm. 

The  quantity  p  has  been  defined  as  the  mass  of  the  swarm  divided  by  the 
cube  of  the  radius.  If  o-  s  cr  (r)  is  the  density  of  the  swarm  as  a  function  of  r  . 
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If  it  is  assumed  that 


(T  e 
o 


then 


2 

P  ={4iro-^/3) 


The  assumed  density  variation  has  been  chosen  conveniently  to  arrive  at  a 
simple  form  for  P(r):  the  density  variation  is  reasonable  since  it  is  an 
exponentially  decaying  function  of  r,  and  of  course  it  is  necessary  that  r' 
be  less  than  3/(2y)  . 

The  differential  equations  now  take  on  the  form 


u"  -  2v'  - (1+u)  =  -  lit  _  p  ^ 

V(l+u)  +  V 

(13) 

2 

v"  +  2u'  -  V  =  -  ,  ■  I  ^  -i-L.  -  p  e  ’ 

I  2  2  o 

VU+u)  +  V 


2  2  2 

where  r  =  u  +  v  and  rr'  =  uu'  +  w'  .  The  Jacobi  integral  can  be  written 


as 


2  p 

(u')^  +  (v’)^  -  (u  V+2u)-  ■  !«  -  (Pq/y)  +  (2  +-£)=  const  , 

V(l+u)^  +  v^  ^ 


(14) 
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and  an  analysis  similar  to  that  for  a  finite  radius  shows  that 

°  l.  (1-r)^ 

The  effect  of  variable  density  is  to  increase  the  density  ratio,  which  is  again 
a  result  to  be  expected  on  physical  grounds. 

(d)  Concluding  Remarks 

The  dominant  result  of  the  simplest  linear  analysis  is  that  the  density 
ratio  must  exceed  3  for  a  stable  swarm  to  exist  in  circular  motion  about  a 
large  attracting  mass.  The  effects  of  finite  size  of  swarm,  eccentricity  of 
orbit,  and  variable  density  all  act  to  increase  the  least  value  of  the  density 
ratio.  For  the  purposes  of  section  II  it  is  reasonable  to  restrict  ourselves  to 
the  linearized  analysis  of  a  spherical  swarm,  but  extend  the  analysis  to  the 
case  of  the  restricted  problem  of  three  bodies.  Again  the  orbit  is  circular 
and  eccentricity  effects  can  be  neglected. 
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II:  MOTION  IN  THE  RESTRICTED  PROBLEM  OF  THREE  BODIES. 

(a)  Equations  of  Motion  In  Rectangular  Coordinates. 

The  restricted  problem  of  three  bodies  will  be  fcarmulated  in  the  usual  way 
[2,  chap.  8];  see  Figure  (3). 


The  plane  of  the  orbit  of  the  two  finite  bodies  m,(x  ,  y  ,  0)  and  mJx  ,  y  0) 

is  the  plane  z  =  0  and  their  mass-center  is  always  at  the  origin  0  of  the 
inertial  system.  The  center  of  the  spherical  swarm  is  at  PCx,  y,  0)  such 
that  PP^  P^  is  an  equilateral  triangle.  Thus  P  is  at  one  of  the  stable 
equilibrium  Lagrangian  configurations  in  the  restricted  problem  of  three  bodies. 
A  general  particle  m  of  the  swarm  is  at  (|,  tj,  &),  near  to  P;  the  maximum 
radius  of  the  swarm  is  small  compared  with  the  orbital  radius  a  ,  It  is 
assumed  of  course  that  the  mass  ratio  j*  of  the  smaller  mass  to  the  total  mass 
is  small  enough  for  P  to  be  a  stable  point. 


#351 


-18 


Simple  geometric  considerations  show  that 
Tcosiip^)  =  (r^  -  t^)/2  , 


(15) 


rsin(<p^)  =  (aN/3)/2  ,  (a  *  r^  +  r^)  , 

and  according  to  the  assumptions  stated  we  also  require 
Xi  =  “(m^/mpx^,  =  -(m2/mj)y2  ; 

=  nM  ,  m^  =  (1  -  |i)M;  (p,  =  m^lA)  ,  (16) 

r^  «  na ,  ^2  *  » 

The  differential  equations  of  motion  of  a  particle  inside  the  swarm  are 


I  a  -Gm^d  -  Xj)/r^  -  Gm^CS  -  -  GP(S  -  x) 

H  a  -Gmj(ti  -  y^}/T^  -  Gm^Cii  -  -  GPCq  -  y)  (17) 

i  a  -Gm^t/r^^  -  Gm^c/r^  -  Gpt 

rf  =  (I  -  x/  +  (,,  -  yj)2  + 

“  (^  ”  *2^^  +  (^  “  )^  +  ^ 


where 


(18) 
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If  however  the  particle  were  external  to  the  nucleus  of  radius  b  the  last 

three  terms  on  the  right  hand  sides  of  equations  (7)  would  be  written,  respectively, 

as 


2  2  2  2 

where  o-  is  the  constant  density  and  r  =  (|  -  x)  +  (ti  -  y)  +  & 

Since  the  mass  near  P  does  not  influence  the  orbital  motion  of  and 
m^  we  have  the  equations  of  motion 

3  3 

x^  =  -  Gm^Cx^  -  Xj)/a  and  y^  *  "Gm^Cy^  -  yj)/a  , 

where  a  =  r^  +  r^,  and  M  =  m^  +  m^  .  Since  x^  =  -(|i/l  -  and 

y^  s  “*((1/1  “  H^)y2  fha  result  is  that  the  differential  equations  can  be  written 

as 

2  3 

where  k  =  GM/a  ,  and  we  assume  that 

x^  =  r^  coskt  =  (1  -  |A)acoskt  , 
y2  =  r2Sinkt  =  (1  -  ^}asinkt  , 

(19) 

Xj  =  -TjCoskt  =  -|i  acoskt. 


yj  =  -Tj  sinkt  =  -pasinkt 
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Since  x  =  rcos(0  +  <p^)  =  rcos(kt  +  <p^)  , 

X  =  r(coskt)  cos(y^)  -  r(sinkt)  sin(<i>^)  , 

(20) 

X  =  {coskt)(r2  -  rp/2  -  (slnkt)(a'N/3/2)  . 

The  results  for  x  and  y  are 

X  =  a  “  ^^coskt  -  slnkt  j  , 

(21) 

y  =  a  sinkt  +  cos  ktj  . 

Finally  we  write 

I  =  a  ^  ~  coskt  -  sinkt  j  , 

Ti  =  a  ^  sinkt  +  coskt  j  ,  (22) 

4  =  aw  , 

which  introduces  a  uniformly  rotating  reference  system;  in  this  new  system  the 
finite  masses  (;i,  1  "  Sf®  fixed  at  the  points  on  the  u-axis  at  distances 
(1  -  p,  K^),  respectively,  from  the  origin.  In  addition  the  coordinates  have 
been  chosen  to  make  u  and  v  small  compared  with  unity  since  (£,  i))  and 
(x,  y)  are  close  to  one  another. 
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(b)  Stability  Criterion . 

The  same  assumptions  will  be  made  here  that  were  made  in  the  linearized 
analysis  of  I,  that  is,  in  the  equations 


u"  -  2v'  -  u  =  (1  -  2|ji)/2  -  -Uii  _  Xu  , 

q 

V  +  2u'  -  V  =  (Vr/2)  -  f  ^  (vt  ,  xv  , 

p  q 


(23) 


w" 


3  3 

-(i-ii)w/p  “tiw/q  -Xw  , 


where  p^  =  (u  +  1/2)^  +  (v  +  ^3/Z)^  +  w^  , 


=  (u  -  1/2)^  +  (v  +  >/372}^  +  w^  , 


(24) 


it  is  assumed  that 


p~^  =  1  -  C3M)  (u  +  vn/I)  , 
q"^  =  1  -  (3/2)  (-U  +  v>/3)  , 


(25) 


to  lowest  order  terms.  If  terms  quadratic  in  u,  v,  and  w  are  neglected  in 
the  expansions  the  equations  can  be  written  as 
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u"  -  2v'  +  (  X  -  3  /4)  -  YV  =  0  , 
V"  +  2u'  +  (  \  -  9/4)  -  YU  =  0  , 


(26) 


w"  +  (1  +  X)w  =  0  , 

where  y  ~  (3^/3/4)(l~2^):  the  last  equation  again  has  only  oscillatory  solutions 
and  is  independent  of  the  first  two  equations.  If  in  the  first  two  equations 
it  is  assumed  that  u  «  Ue''^,  v  s  Ve''^  the  characteristic  equation  is 
determined  by 


(v^  +  X-3/4)  -{2v+y) 

2  ®  0  , 

(2v  -  y)  (v  +  X  -  9/4) 

2 

Again  v  must  be  negative  for  stability,  and  steps  similar  to  that  used  in  I 
lead  to  the  requirement  that 

X  >(3/2)  +  ^  9/16)'  . 

2 

If  ^  s  0,  Y  *  27/16  and  X  >  3,  as  in  I  .  Since  is  small  in  the  earth- 
moon  system  (approximately  1/82)  and  must  be  less  than  0.0385  for  stability 
at  the  points  in  any  case,  it  can  be  assumed  that  is  small  compared 

with  unity.  If  first  order  terms  in  |i  are  retained  the  result  is  that 


X  >  3  -(9/4)^ 
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The  least  value  of  the  density  ratio  is  therefore  changed  by  about  one  per-cent 
in  the  case  of  the  earth-moon  system. 

(c)  Concluding  Remarks 

(i)  The  stability  criterion,  based  upon  a  linearized  analysis,  shows  that 
the  least  density  ratio  of  the  spherical  swarm  is  3  -  (9/4)p.  .  In  the  presence 
of  a  single  body  the  criterion  is  X  >  3  but  it  must  be  recalled  that  the  masses 
used  for  comparison  are  not  the  same  in  both  cases.  In  the  latter  case  we 
have 

>-l  =  <rs/[Mg/(|TTa^)]>3  , 

where  is  the  mass  of  the  earth  and  o-  ^  =  density  of  the  swarm  .  In  the 

former  case  we  have 

\  2  =  cr  g/[(Mg  +  ira^)]  >  3  -  (9/4)  ^  , 

where  is  the  mass  of  the  moon.  Comparison  of  the  two  density  ratios 
shows  that  X  2  =  X  ^  (1  -  n)  , 

or  Xj>3  +  (3/4)|i  , 

in  terms  of  the  mass  of  the  earth.  Since  the  effect  is  small  in  either  case  we 
evaluate  the  density  3M/(4ira^/3)  =  (3M/(4wR^/3)][(R/a)^]  where  M  is  the 
mass  of  the  earth,  R  is  the  radius  of  the  earth  and  a  is  the  lunar  distance. 
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Since  the  first  factor  is  the  mean  density  of  the  earth  (5. 5  gm/cm^}  and  the 
second  factor  (a/R)  is  approximately  60,  the  density  is  approximately 
(1/13,000)  grams/cm  .  This  is  the  least  density  of  a  spherical  swarm  at 
the  lunar  distance,  or  one  located  at  either  L-polnt.  The  density  of  the 
atmosphere  at  sea-level  is  approximately  (1/800)  grams/cm  and  is  decreased 
by  a  factor  of  10  for  each  12  miles  increase  in  altitude,  roughly  [13].  Thus 
the  density  at  12  miles  would  be  (1/8000)  grams/cm^  and  (1/80,000) 

3 

grams/cm  .  The  least  density  of  the  swarm  would  therefore  be  ecjulvalent  to 

the  density  of  the  earth's  atmosphere  in  the  range  12-24  miles  above  sea-level, 

or  well  within  the  re-entry  altitudes. 

(li)  Although  an  explicit  Jacobi  integral  does  not  exist  if  e  #  0  in  the 

restricted  problem  of  three  bodies,  it  is  nevertheless  possible  to  write  as 

follows:  The  same  general  notation  is  adopted  as  in  (Ila),  but  now 
a/1  -  e^l  2 

^2  =  r  =  1  +  ecosv  -  ®  <  1);  v  is  the  true  anomaly,  and 

»  rcosv  ;  =  -  (  (hj:)/(1  -  jt)]cosv  y  r^v  =  h  , 

*  rsinv  ;  y^  *  -  [  (jjj:)/(1  -  |i)  ]  slnv  ;  h^  =  GMa(l  -  e^)(l  -  |x)^  . 

The  true  anomaly  v  is  now  used  as  the  independent  variable  in  place  of  the 
time  t  .  In  the  case  of  a  circular  OTbit  the  differentiations  with  respect  to 
these  quantities  differ  only  by  a  constant  but  in  the  present  case  this  is  not  so. 
Since  r  =  (eh)(sinv)/[a(l  -  e  )]  and  v  =  h/r  ,  we  require  the  change  from  time 
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differentiation  (denoted  by  a  dot)  to  differentiation  with  respect  to  the  true 

anomaly  (denoted  by  a  prime).  If  f  is  a  given  function  we  have  f  s  and 

••  ••  .2 

f  =  f'v  +  f"(v)  ;  with  the  previous  results  this  leads  to 


2 

••  -2eh  sinv  ,, 

2,  3 
a(l-e  )r 


Non-uniformly  rotating  cocnrdinates  are  now  introduced  by  means  of  the 
substitutions 


I  =  Xcosv  -  Ysinv  , 

r\  -  Xsinv  +  Ycosv  , 

C  =  Z  . 

The  finite  masses  are  always  on  the  X-axis,  but  do  not  remain  at  the  same 
locations.  Finally,  since  r  is  now  a  variable  we  introduce  new  distances  so 
that  the  instantaneous  distance  between  m^  and  m^  is  the  unit  of  distance; 
this  changes  with  time  of  course  and  we  write 

Xaw,  Yayr,  Z«2r, 

where  (x,  y,  z)  are  the  new  dimensionless  coordinates.  In  the  subsequent 

2 

transformations  we  use  the  facts  that  r'/r^  *  and  ^  ^  - - — ;  . 

aCl-el  T  f  a(lV) 
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Clearly  if  e  =  0  the  differentiations  with  respect  to  the  true  anomaly  (')  and 
those  with  respect  to  the  time  (*)  are  proportional  to  one  another  and  the 
result  above  reduces  to  the  usual  form  of  the  Jacobi  integral  (2,  p.  281)  • 


#351 


27 


m;  DYNAMICAL  MODEL  FOR  THE  CLOUD-SATELUTES  AT  THE  L .  AND  L, 

4  5 

POINTS  OBSERVED  BY  K.  KORDYLEWSKI  (KRAKOW  OBSERVATORY). 

(a)  Observations  of  K.  Kordvlewskl 

Three  recent  issues  (1961)  of  "Sky  and  Telescope"  have  announced  the 
observations  of  K.  Kordylewski  with  respect  to  cloud-like  satellites  at  the 
L^  and  L^  points,  [14,  15,  16,  17].  In  the  first  of  the  issues  [14]  a  pair  of 
cloud-like  satellites  was  observed  near  the  L^  point,  which  is  the  trailing 
point  60°  behind  the  moon.  The  two  libration  clouds  were  both  near  L^  , 
several  degrees  apart  and  it  was  suggested  that  similar  clouds  could  be  found 
near  L^  ,  The  brief  announcements  in  the  first  issue  were  supplemented  in 
the  August  1961  issue  [15]  where  some  of  the  observational  difficulties  were 
noted.  The  luminous  patches  appeared  to  be  at  least  tivo  degrees  in  diameter; 
two  clouds  about  eight  degrees  apart  were  identified  and  it  was  also  suggested 
that  large  telescopes  be  used  to  locate  individual  meteoroids  in  the  libration 
clouds.  The  third  note,  in  the  December  1961  issue  [16]  confirmed  the  existence 
of  such  clouds  near  the  L^  point  as  well.  The  preliminary  results  indicate 
that  the  L^  object  was  detected  as  a  pair  of  dim  spots,  each  about  five  degrees 
in  diameter  and  nearly  touching.  Some  direct  measurements  were  given  in 
[17]  by  K.  Kordylewski. 

(b)  Basic  Equations  for  the  Dynamical  Model. 

The  dynamical  model  has  been  outlined  in  the  Introduction:  this  consists 

in  the  assvunption  of  a  small  mass  at  or  near  to  one  of  the  L-polnts  (L,  or  L  ) 

^  5 
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and  tenuous  matter  (dust)  in  the  swroundlng  space.  It  will  be  assumed  that 
the  nucleus  as  well  as  the  cloud  near  it  move  in  a  circular  (^bit  about  the 
earth  at  a  distance  equal  to  the  lunar  distance.  The  effect  of  the  moon  will 
be  ignored  in  order  to  simplify  the  presentation.  Since  the  sphere  of  activity 
of  the  moon  [18]  is  too  small  to  enclose  either  L-point  it  is  reasonable  to 
neglect  the  effect  of  the  moon.  More  precise  calculations  show  that  the  effect 
of  the  moon  is  slight  in  this  case,  although  the  moon  does  of  course  determine 
the  location  of  the  L-points  and  the  fact  that  they  are  stable  equilibrium 
points  in  the  earth-moon  system. 

The  small  mass  of  the  nucleus  is  m(x,y)  in  Figure  (4);  a  is  the  lunar 

distance  and  the  plane  of  the  cxbit  of  m 
is  the  (x,  y)  plane.  A  general  particle 
of  the  cloud  is  at  ti,  t)  and  the 
distance  from  the  nucleus  to  any  particle 
is  r^  which  is  small  compared  with  a 
and  r^,  the  distance  to  the  earth  of  mass 
M  at  the  origin;  the  ratio  m/M  =  p  is 
Figure  4  assumed  to  be  much  less  than  unity. 

The  particle  {%,  t|,  is  near  m  but  outside  the  radius  of  nucleus.  The 
general  particle  is  attracted  by  m  and  by  M  but  does  not  influence  their 
motion:  m  is  in  circular  motion  about  the  fixed  mass  M  .  The  equations  of 
motion  can  be  written  as 
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i  =  -GMI/Fj  -  GinCe-x)/r2  , 

Ti  =  -GMti/r^  -  Gintti-y)/r2  , 
i  =  -GM4/rj  -  Gm^/r^  , 

where  G  is  the  gravitational  constant,  and 

rf  -  +  r,^  +  , 

rl  =  (l-x)^  +  (n-y)^  +  , 

X  =  acos(kt),  y  =  asin(kt):  k^  =  GM/a^  , 


-29- 


(27) 


A  rotating  coordinate  system  is  now  introduced  in  which  m  is  a  fixed 
point  on  the  X-axis;  the  unit  of  distance  is  chosen  to  be  a  and  the 
dimensionless  time  t  =  kt  is  introduced  as  the  new  independent  variable, 
with  '  =  d/dr  ; 


I  a  a[X  cos  T  -  Y  sin  t  ]  , 
Tj  *  a[XsinT  -  Y  cos  t]  , 


&  «  aZ  • 


In  addition  we  have  that  m/M  a  p  and 

=  (X-1)"  +  yS  Z"  ,  (p,  Rj,  R^  ,  ^ 


(28) 
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The  final  forms  of  the  differential  equations  are 

X"  -  2Y'  -  X  =  -X/R^  -  P(X-1)/R^  , 

Y”  +  2X'  -  Y  =  -Y/R^  -  pY/R^  ,  (29) 

Z"  =  -Z/Rj^  -  pZ/R^  . 

(c)  Tacobi  Integral;  Discussion  of  Zero  Relative  Velocity  Surfaces. 

If  the  equations  (29)  are  multiplied  by  X',  Y',  Z'  respectively  and 
integrated  the  result  is  the  Jacobi  integral  in  the  form 

s  (X')^  +  (Y')^  +  (Z')^  =  X^  +  Y^  +  -  K'  ,  (30) 

^1  *^2 

where  K*  is  a  positive  constant  in  most  further  discussions.  The  discussion 
of  the  possible  zero  relative  velocity  surfaces  is  simplified  if  cylindrical 
coordinates  are  introduced: 

Z  =  1  'f  pcos  ip  , 

Y  =  psin?»  ,  (31) 

Z  S  z  • 

2 

Terms  up  to  0(p  }  will  be  retained  and  particular  attention  will  be  paid  to  the 
value  z  s  0  and  to  those  values  of  z  that  are  small  compared  with  the  radial 
distance  p  . 
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The  Jacobi  integral  can  be  written  as 

=  3p^  cos^  if,  +  n  -  CK'  +  z^)  ,  (32) 

Vp  +  z 

2 

where  K*  +  z  will  be  denoted  by  K,  a  constant  for  constant  z  .  In 
particular  if  z  =  0  the  Jacobi  integral  can  be  written  as 

=  3p^  cos^  ^  +  (^  -  K  ,  (33) 


which  is  the  form  which  will  be  used  most  frequently.  If  a  particle  is  assumed 
to  be  at  a  given  point  (p^,  ^^)  with  a  given  speed  then  the  constant  K  is 
fixed:  if  V  is  zero  in  Eq.  (33)  the  curves  will  be  the  traces  of  the  zero-V 
siirfaces  (z  s  0),  if  they  exist. 

The  dynamical  equations  of  motion  can  also  be  written  in  terms  of  the 
cylindrical  coordinates.  If  slmplications  are  introduced  which  are  consistent 
with  Eq.  (32)  the  equations  can  be  written  as 


2  2 
p"  “  ~  2py'  =  3pcos  ip  - 


JJL 


,  2^  2.3/2  » 

(p  +  z  )  ' 


»"  +  2pV'  +  2p'  s  -3psin^cos^  , 


(34) 


=  -z 


Jgz. 


,  2^  2,3/2  » 

(p  +  z  )  ' 


2  2  2  2 

where  V  =  (p')  +  (p<i>')  +  (z')  .  However  if  (z/p)  <<  1  the  right  hand 

2  2 

sides  of  Eqs.  (34)  can  be  replaced  by  3pcos  7**  (P/p  );  -3psin^os^;  -z  , 


respectively.  In  the  latter  case  the  last  equation  of  (34)  is  z"  +  z  s  0  which 


32 
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has  only  oscillatory  solutions  and  so  the  motion  of  z  is  known. 

The  dynamical  equations  (34)  of  course  determine  the  actual  motion  of 

each  particle  if  the  initial  conditions  are  prescribed.  However  much  more 

Interest  is  centered  on  the  Jacobi  integral  (33)  which  yields  a  relationship 

between  the  (relative)  position  coordinates  and  the  speed  in  any  realizable 

motion.  That  is,  if  K  >  0  it  is  possible  for  V  to  equal  zero  by  proper  choice 
2 

of  (Pf  ^):  if  the  square  of  the  Initial  speed,  is  sufficiently  large  then 
K  is  negative  and  it  is  not  possible  for  V  to  vanish.  Thus  in  the  latter  case 
no  surfaces  of  zero  relative  velocity  exist  but  in  the  former  case  (K  >  0)  they 
do  exist.  Our  object  is  to  study  such  surfaces,  or  more  precisely  stated, 
the  traces  of  such  surfaces  (z  =  0)  . 

It  may  be  assumed  that  K  has  been  determined  as 

K  =  3j,^cos%  +  (2p/p^) -V^>0  ,  (35) 

and  that  we  have  set  V  >  0  to  obtain  the  equation 

^3  2 

3p  cos  <7-Kp+2ps0  .  (36) 

The  cubic  equation  in  p  has  one  real  negative  root  since  the  leading  coefficient 
and  the  constant  term  are  positive;  this  root  is  of  no  physical  interest.  The 
remaining  two  roots  are  either  conjugate  complex  roots  or  both  positive,  since 
the  sum  of  all  roots  is  zero;  ovir  interest  is  in  the  latter  case,  which  requires 

cos^^  <  K^/(8ip^)  . 


that 


(37) 
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In  the  form  of  the  inequality  (37),  two  cases  of  immediate  interest 

3  2 

arise.  In  the  first,  if  K  >  8ip  ,  then  no  restriction  is  placed  on  <p,  since 

the  inequality  is  always  satisfied;  it  is  assumed  of  course  that  K  must  be 

3  2 

positive  in  all  cases.  In  the  second  case,  if  K  <  8ip  ,  then  there  is  a 
restriction  on  <p  and  not  all  angles  are  possible;  the  admissible  regions  for 


3R^  -  R  +  Y  =  0  , 

where  p  =  R^/i^  and  =  4p^/K^  s  4/9  .  The  cubic  3R^  -  R  +  (2/9)  =  0  has 
the  repeated  root  R  s  -i-  (I/3)  or  p=  ('^)/3  .  The  traces  in  this  case  are 
shown  in  the  sketch  below.  Figure  (6),  and  are  denoted  by  . 


Figure  6 


-34- 


#351 


3  2 

If  K  is  larger  than  8ip  the  traces  are  indicated  by  C^:  for  such  values 

of  K  (p  fixed)  there  is  a  closed  inner  oval  and  a  pair  of  curves  which  extend 

to  infinity  and  reach  the  asymptotic  width  2'V(K/3)  in  all  cases.  It  should 

be  recalled  that  the  unit  of  distance  is  the  lunar  distance  a  .  Thus  as  K 

increases  the  inner  oval  becomes  smaller  in  size  and  the  outer  traces  move 

outward;  the  asymptotic  width  Increases  with  K  .  The  regions  of  positive 
2 

V  are  outside  the  pair  of  curves  and  inside  the  oval.  These  regions  are 

Indicated  for  C.  by  cross-hatched  lines.  The  critical  distance  along 

2  1/3  2/3  2  2 

<»  =  0  occurs  for  C^;  p  =  K/9  =  (3  '  /9)P  '  =  =  p^  .  For  example, 

2  2 

if  a  particle  is  placed  =  0)  at  the  distance  p  with  V  =0,  then  if  p  <  p 

o  c 

2  2 

the  particle  must  remain  Inside  the  oval,  but  if  p  >  p  the  particle  must 

remain  in  the  region  external  to  the  outer  traces  and  cannot  penetrate  to  the 

L-region.  A  particle  placed  at  a  in  Figure  (6)  with  =  0  will  remain  there 

according  to  the  equations  of  motion;  a  particle  placed  at  b  will  move  outward 

but  a  particle  placed  at  c  will  move  inward  (V  =  0}  .  If  the  initial  speed 

—  o 

3 

is  not  zero  then  a  similar  analysis  can  be  made  if  we  require  that  K  be 
2 

at  least  8ip  . 

3 

If  however  is  large  enough  it  is  possible  for  K  to  decrease  and  remain 
3  2 

positive  so  that  K  is  less  than  8ip  but  greater  than  zero.  The  traces  then  change 
in  form  and  only  certain  ranges  of  ^  are  admitted,  as  shown  in  a  typical  case 
in  the  sketch.  Figure  (7). 


2 

The  regions  of  negative  V  are  inside  the  branches  and  the  particles  cannot 
enter  such  spaces.  The  asymptotic  width  of  the  branches  is  still  2v(K/3) 
but  now  K  is  smaller  than  in  the  previous  case  and  so  the  asymptotic  width 
is  much  narrower. 

In  the  present  case  consider  a  particle  in  the  admissible  region  with 

>  0;  see  Figure  (8).  Since  =  3p^cos%  +  ^  -  K  ,  K  «  3p^co8%  +  ^  -V^  • 
o  '^  p^'^o  po 

•'o 

At  point  P,  K  is  fixed  by  the  initial  coordinates  and  speed  (energy).  If  in 
Figure  (8)  the  particle  moves  toward  Q,  V  will  decrease  and  similarly  for 
the  points  R,  S,  ...  which  approach  the  curve  of  zero  relative  velocity. 


Figure  8 
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The  equations  of  motion  ( 34)  of  course  determine  the  actual  path  in  the  relatives 

coordinates.  The  acceleration  is  not  zero  on  the  zero  relative  velocity  curves, 

but  is  normal  to  the  curves  and  directed  outward.  If  some  particles  move  toward 

2 

A,  B,  ...  then  clearly  (Zp/p)  will  increase  and  so  will  V  and  the  particle 
can  move  rapidly  toward  the  L-point  which  is  the  origin  of  the  coordinate  system. 
The  essential  point  is  that  some  particles  may  approach  to  the  V  =  0  curves 
and  remain  there  for  an  appreciable  time  and  a  slow  accumulation  of  matter  can 
take  place  near  such  curves.  If  K  is  smaller  at  P  than  before,  the  curve 
V  s  0  is  thinner  and  inside  the  one  shown  in  the  sketch  above.  Figure  (8).  The 
intersection  point  with  the  axis  is  zp/K  and  this  point  moves  outward.  Thus 
a  narrow  region  of  space  (z  s  0)  acts  as  an  area  where  particles  may  move 
slowly  and  accumulate.  Ultimately  these  particles  will  drift  into  other  regions 
of  space  but  since  there  is  a  steady  drift  of  such  particles  throughout  the  space 
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there  will  always  be  new  particles  streaming  into  the  region  to  replace  those 
which  may  have  moved  away.  This  requires  that  the  particle  speeds  be  close 
to  the  orbital  speeds  at  the  L-points  since  K  is  small  (small  relative  velocity). 
If  the  initial  speeds  are  large  enough  then  no  zero  V  surfaces  can  exist;  this 


is  reasonable  for  it  Indicates  that  the  particles  cannot  be  trapped,  or  brought 

near  to  rest,  in  the  vicinity  of  the  L-point.  This  is  the  case  if 
2  2  2  23 

y>  3p  cos  tp  +  ,  and  therefore  the  right  hand  expression  can  be  considered 

O  O  0 

to  yield  a  minimum  escape  speed  in  the  region  near  the  L-polnts. 

If  z  is  not  zero,  but  still  small  compared  with  then  the  constant  in 

Eq.  (32)  is  merely  increased  slightly.  The  discussion  is  essentially  the  same  as 

2  3  2 

before,  but  now  we  require  that  (K  +  z  )  >  81|3  and  this  also  limits  the  region 
in  the  z-dlrection.  It  is  problematical  if  this  knowledge  could  be  exploited 
observationally  to  determine  an  approximate  value  of  (3;  it  is  not  likely  that 
the  observations  are  clearly  enough  defined  for  this  to  be  the  case,  or  if  the 
form  of  the  clouds  is  in  fact  constant  in  time. 


The  equations  of  motion,  Eqs.  (34),  could  be  used  to  investigate 
further  the  particle  motions  near  the  zero-V  surfaces;  a  simplified  analysis, 
which  is  omitted,  indicates  that  the  particles  will  tend  to  move  away  from  such 
surfaces  (V  =  0)  more  quickly  near  the  L-points  and  less  so  away  from  the 
L-points. 
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(d)  Concluding  Remarks 

If  the  dynamical  model  is  assumed  to  consist  of  a  nucleus  surrounded  by 

dust  particles  moving  near  the  L-regions  then  the  main  tool  in  the  investigation 

of  the  possible  regions  where  such  particles  may  migrate  is  provided  by  the 

Jacobi  integral.  The  constant  K  is  fixed  by  the  initial  conditions.  If  K  is 

negative  then  is  so  large  that  no  zero  V  surfaces  are  possible  and  no 

particle  can  come  to  rest  near  L  (unless  it  strikes  and  is  absorbed  by  the 

3  2 

nucleus).  If  is  sufficiently  small  and  K  exceeds  8ip  then  we  have 

3  2 

the  case  shown  in  Figure  (6),  but  if  K  is  less  than  8ip  we  have  the  case 
shown  in  Figure  (7).  The  asymptotic  width  of  the  curves  is  proportional  to 
•'/iT  and  this  means  that  the  curves  are  much  thinner  in  the  latter  case,  which 
is  the  one  where  the  zero-V  curves  split  into  two  narrow  branches.  Although 
in  the  former  case  a  single  surface  is  possible  this  is  much  wider  and  the 
accumulation  would  not  be  confined  to  successively  narrower  regions.  The 
dynamical  model  assumed  therefore  gives  rise  to  the  possibility  that  faint 
cloud-like  patches  may  appear  in  the  region  of  the  L-points. 
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APPENDIX:  STUDY  OF  THE  STABIUTY  AND  MOnON  OF  THE  TWO-PARTICLE 
MODEL. 

In  the  Introduction  the  circular  motion  of  two  equal  infinitesimal  particles 

of  mass  m  was  considered  in  relation  to  a  fixed  mass  M  .  The  value  of 
3  3 

X  =  mb  /(4Ma  }  is  3  if  (a/b)  is  much  less  than  unity.  If,  for  example, 
m  =  2500kg,  then  2m  =  5000kg  or  about  11,000  pounds.  If  two  such  masses 
are  in  close  circular  orbit  about  the  earth  it  is  interesting  to  compute  how 
close  two  such  bodies  must  be  for  p  to  equal  3  .  If  the  distance  b  is  kR  , 
where  k  is  a  constant  multiplier  and  R  is  the  radius  of  the  earth,  and  use  is 
made  of  the  fact  that  the  mean  density  of  the  earth  is  5.5  gm/cm  then  the 
distance  between  the  two  particles  (2a)  is  given  by  2a  s  41. 5k  in  centimeters 
for  the  given  data  (m  s  2500kg}  •  Thus  at  a  distance  of  10  earth  radii 
(40,000  miles),  2a  =  415cm  or  4.15  meters  or  approximately  13  feet  . 

These  figures  would  change  if  the  distances  and  densities  would  be  those  for  the 
moon,  since  the  density  is  lower  and  the  radius  is  smaller. 

The  question  of  the  stability  of  such  an  assumed  motion  can  be  investigated 
in  a  simplified  way  if  we  restrict  ourselves  to  the  orbital  plane  of  motion; 
see  Figure  (9). 

(L  >0)  , 

»  (02  “  ®i)  • 


Figure  9 
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The  equations  of  motion  in  polar  coordinates  are 


•  2  2  3 

"  “GM/fj  -  (Gm/L  )(rj  -  r^cos?*)  , 


ri0j  + 2^101  =  (Gm/L  )(r2Sin<i>)  ; 

~  ^2^2  ~  ~  (Gm/L^)Cr2  -  r^cos^)  , 

’^2®2^^^2®2  =  (Gm/L^)(rjSin  . 


(38) 


In  the  case  of  uniform  circular  motion  it  is  assumed  that 


r^  =  Rj  =  constant, 


r^  =  R2  *  constant  , 


e  =  0-  =Cwt)  , 


u  =  constant 


(39) 


which  implies  that  ?»  =  0  and  L  »  Il^-Rjl  =  d  >  0;  the  second  and  fourth 
equations  in  Eqs.  (3^  are  identically  satisfied  and  eliminaUon  of  between 
the  two  remaining  equations  leads  to 

.3,.2 


^  (RjR2)^(Ri  +  R2) 


(40) 


which  is  equivalent  to  the  relation  (1)  .  The  following  two  expressions  for 
will  be  of  use  in  subsequent  wcric: 
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=  (GM/Rj^  -  (Gm/d^)(R2-Rj)/Rj  ,  (41) 


«^  =  {GM/r|)  +  CGm/d^)(R2-Rj)/R2  .  (42) 


In  the  slightly  disturbed  motion  it  is  assumed  that  each  of  the  4  coordinates 
may  be  varied;  that  is , 


*^1  "  ^  ^2^  ' 
I  ^^2  ~  ^2^^  ^  ^2^  * 


0j  s  (a>t  4  , 

®2  =“t+“2  » 


(43) 


where  pp  p^,  are  dimensionless  quantities.  It  is  also  assumed  that  tp 

is  so  small  that  cos  s  1  and  that  sin  ^ ’s  q,  a  •  a^;  only  terms  linear  in 

p,  a  will  be  retained  in  the  stability  analysis.  Finally  we  set  a>t  =  t  with  ' 

denoting  d/dr,  and  approximate  the  binomial  expansions  in  the  usual  way: 

-2  -2 

r^  s  R^  (1  -  2p^)  and  similarly  for  subscript  Z,  and 


The  first  dynamical  equation  can  be  written  in  the  following  way,  after 
simplifications  have  been  made: 


Pi 


'1  J 


+  2p, 


QM.Qm 


(44) 


■44- 


#351 


2 

The  relation  for  u  ,  Eq.  (44  permits  this  result  to  be  simplified  to  the 
following  form  if  Gm/(«^d^)  =  6^  and  (Rj/R2)  ®  ®2  ’ 


p"-3pj-2aj  =  -2{6j/62){p2-P2)  • 


(45) 


The  second  dynamical  equation  can  easily  be  simplified  to  the  form 


a'>  +  2p|  =(6^/62) («2“«i)  •  (46) 

The  remaining  pair  of  dynamical  equations  can  be  simplified  in  the  same  way 

2 

if  the  alternative  expression  for  <<>  ,  Eq.  (42),  is  used.  The  four  equations 
can  be  written  as 


pj  ”  3pj  -  2arJ  *  “2(8j/62)(p2  “  Pj)  *-3(p2  "  f^) 


«;•  +  2pi  =  (6i/62n«2  "  “1)  *  (3/2)(«2  -  «j)^ 
P2  “  ^P2  -  2“2  =  ®2>tp2  ’  Pi?  =  ^<P2  “  Pl> 


(47) 


"2  ^  2P2  *  ^®1  ®2^t*2  "  "l^  '  (3/2)(«2  -  “1) 

where  the  second  set  of  right  hand  sides  are  approximations  to  those  which 
precede  them  since  (R2/RJ)  =  ^2  close  to  unity  and  (Gm)/(w^d^)  s  6j 

is  nearly  3/2  under  the  same  assumptions.  If  the  simplified  forms  are  taken 
the  above  set  of  four  equations  can  be  decomposed  into  the  following  two  pairs 
of  equations  by  straightforward  additions  and  subtractions.  The  first  pair  of 
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(P2+  Pi)"  "  3(p2  +  pj)  -2(0-2  +aj)'  =  0 
(«2  +  orp"  +  2(p2  +  Pj)'  =0 

(P2-Pir-9(P2-Pi)-2(«2-«i)'=0’ 

(*2  -  flfj)"  +  3(0-2  "  “P  ^**2  " 


•45 


(48) 


(49) 


equations  relates  to  the  motion  of  the  mass-K:enter  of  the  pair  of  particles  and 
the  second  pair  to  the  relative  motion  of  the  particles.  To  determine  the 
solutions  let 


UT  UlT 

Pi  +  P2  *  Ae  ,  or,  +  a,  s  Be  , 


1  ■  2 


and 


«  P-'*’  ^  P"’’ 

Pi  -  =  Ce  ,  Oi  -  a2  =  De  , 


in  the  two  sets  of  equations  respectively.  The  two  characteristic  equations 
are 


(P^  -  3)  -Z\L 

(v^  -  9)  -2v 

s  0  and 

y  2 

2p  P 

2v  (v^  +  3) 

0  , 


P  =0,  -1 


or 


and 


=  1  *  's/28 


(50) 
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Thus  the  disturbed  motion  is  unstable  both  with  respect  to  the  motion  of 
the  mass-center  and  the  relative  motion  with  respect  to  the  mass-center.  The 
complete  solution  for  all  four  quantities  can  now  easily  be  carried  out,  but 
more  direct  interest  is  in  the  solution  for  the  relative  motion;  this  is  the 
motion  which  characterizes  the  separation  of  the  particles.  Here  we  have  the 
relations  C3(v^  -  9)  =  (2v)D  for  the  coefficients.  If  s  1  +  n/28  ,  v  =  *2. 51 

and  C/D  »  +  1.86;  if  =  1  -  ^fz8  ,  v  =  A  i2.08  and  C/D  =  +  iO.156.  If 
new  coefficients  are  introduced  as 

Dj=(Ej+E2)/2,  D^  =(Ej-E2)/2  and  Dj  =(E3+1E^)/(21)  ,  D^  =  (E3-iE4)/(2i)  , 
the  solutions  for  the  motion  of  separation  can  be  written  as 

a2  ■  «i  =  El  cosh(2,  51t)  +  E^  sinh(2.5lT}  +  E3COs(2,08t)  +  E^sin(2,  08t),  (51) 

p^-  Pj  =  ”  1. 86EjSinh(2.  51t)  - 1,  86E2COsh(2.  51t)  -  0. 156E3Sin(2.  08t)  -  0, 156E^cos(2.  08t)  , 

(51) 

and  the  solutions  for  the  motion  of  the  mass-center  can  be  written  as 


«2  +  “i  =  (^1  +  +  F3COs(t)  +  F^sin(T)  , 

P2+Pj  =  (2/3)F2  +  (l/2)F3Sin(T)-(l/2)F^cos(T)  , 


(52) 


in  terms  of  new  constants;  there  are  8  arbitrary  constants  since  the  position 
and  velocity  of  each  particle  can  be  prescribed.  At  t  s  0  all  4  coordinates 
may  be  assumed  to  be  equal  to  zero,  that  is,  or.  *  a  =  0  and  p,  p,  s  o  and 

Z  1 
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only  4  constants  remain  in  the  solutions: 

a  -  a  =  E,(cosh  2,  51  t  -  cos  2.  08t)  +  E_  (sinh  2.  51t  -  11. 9  sin  2.  08t)  , 

2  11  ^ 

=  Ej(-1. 86sln  2.  51t  +  0. 156sin  2.  08t)  +  0. 186E2(-cosh  2,51r  +  cos  2.  08t)  , 
0-2  +  «i  =  Fj(l  -  cos  t)  +  F^Ct  -  (4/3)  sinr)  ,  (53) 

P2  +  Pj  =  -(2/3)F2(1  -  cost)  -(1/2)Fj  sin  t  . 

The  equations  may  be  solved  under  various  initial  conditions. 

For  example  we  may  set  (p2  -  Pj)'  =  0  and  -  a^)'  s  and  determine 
the  time  t  =  t’X  at  which  (a,  -  a  )  vanishes.  This  requirement  leads  to  the 
equations 


0.044q  [sinh  2,51t  -  11.9  sin  2,08t]  ,  (54) 

2  1  O 

(54) 

P2  -  Pi  =  0. 186(0.  044q^)  [-cosh  2,  51t  +  cos  2.08t]  , 

and  if  a,  *  a,  at  t  =  t*  ,  (  p-  -  p,)  *  P^,  a  given  number,  then 
2  1  ''•^2  10' 

sinh  2.51t*  »  11.9  sin  2.  08t*  . 

The  value  of  t*  is  uniquely  determined  from  this  equation  and  corresponds 

to  the  condition  second  equation  of  Eq.  (54)  fixes  the 

value  of  q  for  this  motion  to  take  place.  That  is,  q  is  now  fixed  so  that 
o  o 
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at  T  =  T*  the  angle  tp  will  be  zero  and  the  distance  of  separation  will  be 
zero;  in  order  for  this  to  hold  under  the  assumed  initial  conditions,  only 
one  value  of  is  possible. 


